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Abstract 

The use of an external transverse magnetic field to trigger and to control electron self-injection 
in laser- and particle-beam driven wakefield accelerators is examined analytically and through full- 
scale particle-in-cell simulations. A magnetic field can relax the injection threshold and can be 
Qh used to control main output beam features such as charge, energy, and transverse dynamics in the 

{A 

O ion channel associated with the plasma blowout. It is shown that this mechanism could be studied 

•i-H 

^ using state-of-the-art magnetic fields in next generation plasma accelerator experiments. 
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Plasma based acceleration (PBA) has the potential to lead to a future generation of com- 
pact particle accelerators [IH3] for high energy physics applications, and also to compact light 
sources for medical applications [4] and material science. In state-of-the-art laser wakefield 
acceleration (LWFA) experiments, relativistic plasma waves capable to self-inject and accel- 
erate 1 GeV class electrons [5] are excited in the so-called bubble or blowout regime [6l [7j . 
In contrast, the self-injection conditions are not easily met in current plasma wakefield ac- 
celeration (PWFA) experiments, since the bubble radius is not as large. Nevertheless, in 
recent PWFA experiments, some electrons had their energy doubled from 42 GeV to 85 GeV 
in 85 cm [8j. 

The output controllability is a major challenge for the use of PBA in several applications. 
To this end, new concepts have emerged such as counter- and cross-propagating laser pulses 
for LWFA [9], short plasma down-ramps [10], ionization induced trapping [11], and evolving 
bubbles [12]. With these methods, the charge and energy of self-injected bunches can be 
adjusted. 

In this Letter we propose a novel scheme that uses static transverse magnetic fields 
to trigger and to control the self-injection in the LWFA or PWFA. The trapping occurs 
off-axis leading to synchronized betatron oscillations, which could improve the quality of 
x-ray emission by the magnetically injected electrons [13]. The output energies can also be 
controlled by adjusting the longitudinal injection position. The scaling law for the magnetic 
field induced injection is determined with the appropriate Hamiltonian, and illustrated with 
particle-in-cell (PIC) simulations in OSIRIS [2]. For the next generation PBAs aiming at 
producing multi-10 GeV electron bunches in controlled injection scenarios [2|3], our scheme 
requires external B-fields as low as 5 T. 

The use of magnetic fields in PBAs was first explored in the so-called surfatron [T5] . 
The role of magnetic fields in PBAs was further examined in [16], while Ref. [17] showed 
that external longitudinal magnetic fields increase the trapped charge. Our scheme uses an 
external magnetic field perpendicular to the driver velocity to generate controlled off-axis 
injection bursts. 

Self-injection can be investigated through the Hamiltonean H = ^m 2 e c A + (P + eA/c) 2 — 
ecj) for background electrons, where P = p — eA and p are the canonical and linear momenta, 
— e the charge of the electron, and A and (j) the plasma vector and scalar potentials. Unless 
stated, normalized units are adopted henceforth. Durations are normalized to the inverse of 



the plasma frequency uj p = ^/47m e 2 /m e , lengths to c/uj P) velocities to c, momenta to m e c, 
and charge to e, with n the background plasma density, and m e the electron mass. Vector 
and scalar potentials are normalized to e/m e c 2 and to e/m e c. The normalized B-field is 
given by u) c /^ P = e|B|/m e cc;p, and uj c is the cyclotron frequency. 

In the co-moving frame (x — x,y — y,£ — v^t — z,s = z\ the Hamiltonian is % = 
H — Vcf)P\\, where v^ is the plasma wave phase velocity, and Py is the longitudinal canonical 
momentum. In this frame dl-L/dt; = (1 — v z /v^)~ (v • dA/ds — d(j)/ds) [HJ [18], where v z 
is the electron velocity parallel to the driver velocity. Thus, defining the wake potential 
as ip = (j) — v$A z [IB], and since v z = v^ for trapped electrons, 1 + At/j = 7/7? — Al-L 
at the instant of injection, where 7^ = (1 — ^|) -1 ^ 2 is the relativistic factor of the bubble 
phase velocity, and Aip = ijjf — ip i (AH = Hf — Hi) is the difference between the ip (Ti) 
of the electron at its trapping and initial position. This trapping condition is consistent 
with the conditions from other electron self-injection mechanisms [HJ H21 HEJ EO] , and is also 
valid when external fields are present. Considering that ip ey:t is related to the external field, 
and denoting the plasma wave contribution to ip by ^ pl , the trapping condition becomes 
1 + A^ 1 = 7/72 -Al-L- A^\ 

The trapping condition derived above shows that self-injection can be relaxed by the 
presence of external fields. In the absence of external fields (Aip ey:t = 0) and in the quasi- 
static approximation [20J, such that AH = 0, self-injection occurs when A?/; pl approaches 
A^ pl = —1 [HI [20j EI]. This condition can be relaxed (i.e. lower |A^| can still lead to 
injection) during the expansion of the wakefield [12] which contributes with finite AH. In 
the presence of external fields, self-injection may be controlled by acting on Al-L and Atp eyit . 
To illustrate the latter mechanism, an external constant magnetic field B y pointing in the 
positive y-direction is considered, and described as ip e ^ = — A e ^ = B y x. We also assume 
that the B-field rises from zero to B y0 in a length L ramp , is constant for L flat , and vanishes 
in L ramp . 

We start by investigating trapping in the uniform B-field region, for which the trapping 
condition reduces to 1 + A^ pl = —B y0 Ax when 7^ —} 00, and where Ax = Xf — Xi is the 
difference between the final and initial electron position in the x direction. A conservative 
threshold B field for injection can be retrieved assuming that A?/; pl = 0, and that the wake 
is relatively unperturbed by the B-field (i.e. uj c /uj p «1). In LWFA matched propagation 
regimes, the B-injected electrons originate at a distance Xi ~ r^ from the axis, and they 



are trapped with Xf ~ [21j. Thus, for a trapped electron Ax ~ —r^sinO = — 2y / a^sin#, 
where r^ = 2^/a^ is the blowout radius [7], # the angle between the plane of the trajectory 
of the electron with the B-field, and a is the normalized laser vector potential. Then, the 
threshold B-field for injection becomes: 

B t = - • (1) 

2sm9^ v ; 



Eq. (fTl) illustrates that the external B-field leads to localized off- axis injection in a well 
defined angular region. We note, however, that since typically AiJj < — 1, the assumption 
Ai/j = significantly overestimates the required B-field for injection. In fact as we will 
show in 3D OSIRIS PIC simulations of the PWFA, the threshold B-field for injection is 
significantly lower and can be within the reach of current technology [19] . 

The localized injection in the transverse x-y space can also be interpreted in terms of 
the Larmor rotation of the plasma electrons. Since backward moving electrons in the z- 
direction rotate anticlockwise for B y0 > 0, they are bent towards the axis for x > 0, entering 
the accelerating and focusing region of the bubble with larger p\\ ) thus facilitating self- 
injection. For x < 0, electrons move away from the axis, preventing injection. Furthermore, 
the electrons moving in the plane x ~ are much less disturbed. Hence, the B-field leads 
to off-axis injection in a defined angular region. 

If B y 3> B t the plasma wave structure can be significantly modified, which may reduce 
|A^| and |Ax|, and suppress injection. In such cases, however, injection can still occur 
because the B-field depends on £ (contribution of AH, and A^ ext to the trapping condition). 
Specifically, as the B-field begins to decrease in the down ramp region then ip pl returns 
to its unperturbed state, and r^ gets larger, effectively lowering v^. In the down ramp B- 
field regions injection also occurs off- axis and in an angular region that is identical to that 
associated with the uniform B-field regions. We stress that for the scenarios presented here, 
the majority of the charge is injected in the B-field down ramps. This is the most important 
mechanism in the present configuration, although for sufficiently large L ramp the mechanisms 
associated with constant B-fields may also lead to self-injection. 

According to the trapping condition in the presence of external fields, higher injection 
rates occur for higher plasma wave expansion rates, or equivalently for higher B-field down- 
ramp gradients, i.e. for shorter ramps or higher B's which provide a larger \v-dA/ds— d<j)/ds\ 
and hence larger \AT~L\. Physically, both tend to lower the phase velocity due to the accordion 



effect [22]. We note that L ramp can change by modifying the spatial extent of the B-field, 
and by varying n . 

In order to illustrate this controlled injection mechanism we have performed a set of 
3D PIC simulations in OSIRIS [Hj. We first consider the LWFA scenario, illustrated in 
Fig. nl The simulation box with dimensions of 24 x 24 x 12 (c/uj p ) 3 moves at the speed 
of light, and is divided into 480 x 480 x 1200 cells with 2x1x1 particles per cell. The 
ions form an immobile neutralizing background fluid. A linearly polarized laser pulse with 
central frequency ujo/uj p = 20 was used, with a peak vector potential of a = 3, a duration 
^p^fwhm = 2^/a^, and a transverse spot size (Wo) matched to the pulse duration such that 
Wo = ctfwhm The plasma density is of the form n = rio(z) (1 + Anr 2 ) for r < 10 c/uj p 
and n = for r > 10 c/uj p with An = 4/W^ being the linear guiding condition, and where 
rio(z) is a linear function of z which increases from n = to n = 1 in 50 c/u) p ensuring 
a smooth vacuum-plasma transition. At the point where the plasma density reaches its 
maximum value, a constant external field pointing in the positive y-direction rises with 
B^ = uj c /uj p = 0.6sin 2 [7r(z - z )/(2L ramp )], with L ramp = 10 c/u p and z = 62 c/u p , 
it is constant and equal to Bff = 0.6 for L flat = 40 c/u) p and drops back to zero with 
B^ = 0.6sin 2 [^(z - Zi)/(2L ramp )] with z 1 = L ramp + L flat + z . For a laser with central 
wavelength of 800 nm, these parameters correspond to a 0.473 J laser pulse, with Wo = 
8.68 /xm, and t FW hm = 28.9 fs. Moreover, n = 4.5 x 10 18 cm" 3 , and B y0 = 407 T. 
While this simulation clearly identified key physical mechanisms of the B-field injection using 
reasonable computational requirements, additional 2D simulations (not shown) revealed that 
self-injection assisted by magnetic fields can be achieved in the LWFA within state-of-the-art 
magnetic field generation technology [T9j . 

Fig. [IJi shows the electron density, the laser projections, and electrons above 10 MeV 
(spheres). Each quantity is plotted when the front of the laser has just left the B-field, 
at uj p t = 120, when the laser evolution [23] can be neglected, corresponding to the wake 
lying entirely within the downramp region. Fig. [l)>d shows how the B-field can lead to 
self-injection. Although the B-field can still decrease r\> in the y-£, x=0 plane, the wake 
remains symmetric. This can be seen as follows: when an electron is expelled sideways the 
B-field provides no extra force. However, as the B-field within the wake bends the trajectory 
backwards in z, the external B-field (in y) produces a force in the negative x direction. This 
motion then provides a force in the positive z thereby decreasing r&. On the other hand, the 
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FIG. 1. 3D OSIRIS simulation of a magnetized LWFA. a. Electron density iso-surfaces in the 
uniform B field region, b-d. Density profile associated with the y = 0, £ = 0, and x = planes of 
a revealing the off-axis injection. Self-injected electrons (darker dots) are closer to the bubble axis, 
while the electrons farther from the bubble axis escape from the trapping region (lighter dots). 
e. Phase-space of the plasma electrons and spectrum (solid line), showing a quasi- monoenergetic 
(~ 6% FWHM spread) electron bunch, f. Trajectories of the self-injected beam particles, with the 
propagation axis shown by the dashed line. The B-field profile is also represented. 



wake is asymmetrically modified in the x-z (y=0) plane. Electrons moving backwards in z 
feel an downward force from the external B-field. This external force reinforces (reduces) the 
focusing force for electrons with x > (x < 0). This leads to the sheath structure seen in 
Fig. [pL As predicted from the trapping condition, electrons with Aip ext < or equivalently 



Ax < are more easily trapped, thus guaranteeing off axis injection. 

The localized trapping provided by the external B-field is also seen in Fig. [IJi-c where all 
electrons above 10 MeV are shown as spheres (dots) at uj p t = 120. The electrons all reside 
with x > and are localized in y as well. Electrons which eventually reside outsize the 
wake (large x) are actually defocused by the external B-field. These energetic electrons now 
move forward in z such that the external B-field leads to a force in the positive x direction. 
Electrons outside the wake do not feel the focusing force of the wake and therefore are lost. 

In the B-field down-ramp, r^ increases, effectively decreasing v^ and leading to an ad- 
ditional and stronger injection. The energy spectrum of the electrons at later propagation 
distances are shown in Fig. [J, revealing that a quasi mononergetic beam is formed at 
ujpt = 270. The total charge of electrons above 50 MeV is 0.15 nC. The streaks in Fig. fife 
correspond to the trajectories in Energy-^ phase space, and are plotted during the time 
between 240-270 uj~ 1 . The color corresponds to their initial value of £ (£i n j), then show- 
ing that most electrons originate in the downramp region 112 < £ mj < 122. Fig. [jj shows 
that self-injected electrons start with x > forming a beam that has a centroid executing 
betatron oscillations. 

In Fig. [2] this scheme is applied to PWFA. The energy of the electron beam driver is 
30 GeV, and its density profile given by n& = n^exp (— x^_/(2<r[)) exp (— ^ 2 /(2a^)), where 
a± = 0.3 c/uj P) a z = 0.5 c/uj P) and n^ = 8.89 n . For n = 10 15 cm -3 , this corresponds to 
a± = 50.4 /im, a z = 84 /xm and to a total number of 3 x 10 10 electrons, close to the SLAC 
electron beam. The simulation box is 12 x 12 x 16 (c/uj p ) 3 , divided into 480 x 480 x 640 
cells with 2x1x1 particles per cell for the electron beam and background plasma. The 
B-field profile is similar to the LWFA case with B^ = 0.55 oj c /u p . For n = 10 15 cm -3 , 
this corresponds to L flat = 6.8 mm, 27 amp = 1.7 mm, and £>Sf = 5.5 T, within current 
technological reach [19]. 

Injection is absent in the unmagnetized PWFA (Fig. [2^i-b). In the magnetized case, we 
observed that self-injection occurs only in the B-field downramp, leading to the generation 
of a 13 pC electron bunch with 0.1 c/uj p = 16.8 /im long, and 0.2 c/uj p = 33.6 fim wide for 
no = 10 15 cm -3 . As in the LWFA, the deformed structure of the wave only traps plasma 
electrons located in a narrow angular region (Fig. [2]>d), resulting in synchronized betatron 
trajectories of the B-injected electrons. We note that the external field also deflects the 
driver by an angle 5 ~ B^L/^ where L is the total B-field length, and 75 is the electron 



beam relativistic factor. In our scenario, as 5 ~ 6 x 10 -4 , the driving beam deflection is 
negligible, but can still be corrected by placing additional identical external B-fields pointing 
in the opposite direction. In fact, separating the additional fields by the betatron wavelength 
can further enhance synchronized injection. 
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FIG. 2. 3D OSIRIS simulation of B-injection in PWFA (2D slices represented), a. Electron 
density, and b. transverse density slice at the back of the bubble in the unmagnetized case. c. 
Electron density after the B-field down-ramp and d. Transverse density slice at the back of the 
bubble, revealing the trapped particles. 



Resorting to 2D simulations, the role of the B-field amplitude, and down ramp length was 
examined for the LWFA case (Fig. [3]), starting with the parameters associated with Fig. [TJ 
but using uJo/ujp = 50. Fig. [3] supports our predictions, showing that the self-injected charge 
decreases with L ramp (Fig. |3k), and increases with Bff (Fig. |3b). The simulations also 
showed that the self-injected beam radius, energy spread, and beam emittance generally 
lowers with the applied B-field amplitudes and for shorter ramps. Moreover the beam 
duration increase for larger B-fields, and for larger field downramps. The injection shutdown 
occurs for L ramp larger than a few hundred electron skin-depths, corresponding to the cm- 



mm scale for plasmas with no = 10 — 10 cm . Finally, these simulations indicate that 
the threshold B field for injection in the down-ramps is similar to that associated with the 



uniform B-field region, corresponding to £>^ xt = 1 — 100 T with n = 10 
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FIG. 3. Parameter scan for the injected charge in the first bucket retrieved from 2D simulations, a. 
Down-ramp length scan, revealing that charge decreases with L ramp b. Peak B-field scan, revealing 
the charge increases with £>^ xt . 

In conclusion, we explored a novel controlled injection mechanism valid for both LWFA 
and PWFA. This scheme has the potential to generate high-quality beams, since it allows 
for the tailoring of the injection time period and azimuthal range. It is possible to generate 
0.1-1 nC-class electron bunches depending on the plasma and beam density parameters. The 
plasma, and B-field dependence for the injected charge, and self-injection were estimated 
analytically and are consistent with the results from PIC simulations in OSIRIS. 
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